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Abstract. We investigate the solutions to the Lorentz-Dirac equation and 
show that its solution flow has a structure identical to the one of renor- 
malization group flows in critical phenomena. The physical solutions of the 
Lorentz-Dirac equation lie on the critical surface. The critical surface is 
repelling, i.e. any slight deviation from it is amplified and as a result the so- 
lution runs away to infinity. On the other hand, Dirac's asymptotic condition 
(acceleration vanishes for long times) forces the solution to be on the critical 
manifold. The critical surface can be determined perturbatively. Thereby 
one obtains an effective second order equation, which we apply to various 
cases, in particular to the motion of an electron in a Penning trap. 



The Lorentz-Dirac equation governs the motion of a classical charge in 
prescribed external electromagnetic fields and including radiation reaction, 
i.e. the loss of energy due to radiation. In standard relativistic notation it 
reads [|1| 

mv" = eF"" {x)v^ + (e Vevrc^) [v^" - \ v^vxv^"] . (1) 

Here m is the (experimental) rest mass of the particle with charge e. x'^{s) is 
the world line and v'^{s) = x'^{s) the velocity of the charge parametrized in 
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its eigentime s. The particle is subject to time-independent external fields 
as given through the electromagnetic field tensor F'^'^. The first term in (1) 
is the Lorentz force while the second term describes the radiation reaction. 

One obvious issue is to understand how the Lorentz-Dirac equation is 
related to the Maxwell-Lorentz equations with a suitable ultraviolet cut- 
off. This problem was studied extensively by Abraham, Lorentz, and many 
others, cf. || for a detailed account. In his famous paper Dirac 
circumvented the issue through a somewhat delicate splitting of the fields 
generated by a point charge. The, to our knowledge most complete formal 
derivation of (|1|) has been worked out by Nodvik [Q. Some rigorous results 
are |0, ||, 0]. For the purpose of this letter we regard the Lorentz-Dirac 
equation as given. 

As noted already by Dirac, Eq. (|l]) has runaway solutions which grow ex- 
ponentially in time, simply because for F'^'^ = and in the approximation of 
small velocities we have mv = (e^/67rc^)v. Dirac 0, reemphasized by Haag 
P), postulated that the physical solutions to (|l|) must satisfy the asymptotic 
condition lim v'^{s) = 0, which, as extra bonus, is a substitute for the miss- 

s— >oo 

ing initial condition x(0). The validity of the asymptotic condition has been 
tested only in explicit cases |^, |TU[. With a general external field tensor 
Ff^^ the solution behavior of (|l]) might be complicated and should expected 
to be chaotic. Physical and unphysical solutions might be thoroughly mixed. 
Thus in principle, for given x(0), x(0), there could be many solutions satisfy- 
ing the asymptotic condition. Which one to pick then? On a more practical 
level, one would like to have a reliable numerical scheme not hampered by 
the instability of physical solutions. 

The purpose of this letter is to explain that the solution fiow of the 
Lorentz-Dirac equation has a structure familiar from the renormalization 
group fiows in critical phenomena. The physical solutions lie on the critical 
surface, which contains attractive fixed points whose location depends on 
F'^". Slightly off the critical surface, the solution grows exponentially fast, 
so to speak it ffows to the high, resp. low, temperature fixed point. Our 
observation has two important implications. (1) The critical manifold is 
actually a surface of the form x = /i(x, x). Thus, for given initial conditions 
x(0),x(0), there is exactly one solution on the critical surface and, as to 
be shown, it satisfies the asymptotic condition. (2) There is an effective 
second order equation, given below, which governs the motion on the critical 
surface. Thus the initial value problem is restored and the equation can 
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esasily be solved numerically. We will demonstrate the predictive power of 
the second order equation by a few examples, still handled without numerical 
integration, the physically most relevant of which is the motion of an electron 
in a Penning trap fl^. (This system was pointed out to us by Wolfgang 
Schleich). 

In all applications the radiation reaction is a small correction to the 
Lorentz force equation, which means that the radiation reaction term, the 
highest derivative in (|ID, carries a small prefactor. Differential equations of 
such a type have been studied extensively through singular (or geometric) 
perturbation theory |]13|, |T^ , which is closely connected to the theory of cen- 
ter manifolds. The application to the Lorentz-Dirac equation is a little bit 
messy and has been carried out in . Rather than trying to summarize these 
results, we believe it to be more instructive to illustrate the basic features of 
the method by using a fictituous mathematical example. 

Let us consider then the ordinary differential equation of the form 

X = gix, y), ey = y- h{x) , (2) 

x{t) G M, y{t) G M. We want to understand the behavior of solutions for 
small e. If we simply set e = 0, the second equation reduces to y = h{x) 
and therefore x = g{x, h{x)). The ambient phase space has disappeared and 
the motion takes place only on the one-dimensional surface {y = h{x)}. On 
the other hand we can go over to the slow time scale t,t = s^^t. Denoting 
differentiation with respect to r by ', (|^) reads 

x' = eg{x, y), y' = y - h{x) . (3) 

Setting now e = 0, yields x' = 0, i.e. x(r) = xq and y' = y — h{xo). 
Thus the surface {y = h{x)} consists exclusively of repelling fixed points. 
If 1/(0) 7^ h{xo), the solution grows exponentially. In this sense the surface 
{y = h{x)} is critical. The main result of geometric singular perturbation 
theory is that for small e the critical surface persists and is of the form 
{y = he{x)}. On the critical surface the motion is governed by 

X = g{x,h^{x)) . (4) 

If x(0), y{0) are off the critical surface, the solution to (Q) diverges exponen- 
tially with rate 1/e. 

Of course, abstractly only the existence of h^ is asserted. Its concrete 
form must be extracted from (^. Fortunately we are allowed to determine 
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hi; perturbatively (which is not the case for individual solutions). We make 
the ansatz y = h^{x) = hoi^x) +£hi{x) + Oi^e^) and insert in (||). This results 
in eij = ehQ{x)x + 0{e'^) = eh'(^{x)g{x, hQ{x)) + 0{e'^) = ho{x)+ehi{x) — h{x) + 
0{e^). Comparing orders of e yields ho{x) = h{x), hi{x) = hQ{x)g{x, ho{x)). 
Inserting in (^), we have the effective equation of motion 

X = g{x, h{x)) + edyg{x, h{x))h'{x)g{x, h{x)) , (5) 

valid up to an error of order e'^. 

Returning to the Lorentz-Dirac equation, by the same argument it has a 
repelling critical surface of the form {x = /i(x, x)}. To understand the motion 
on the critical surface we fix one inertial frame and denote the position and 
velocity three-vectors by r(t), u(t). If is the external electrostatic potential, 
then the energy balance reads 

— [ymc^ + e0(r) — (e^/67rc^)7^(u ■ u)] 

= -(eV67rc3)[7'u' + c'Vlu ■ u)'] (6) 

with 7 = 1/ a/1 — u^/c^. We integrate both sides of in time. Since on 
the critical manifold the Schott term — (e^/67rc'^)7^(u ■ ii) is bounded, we 
conclude that 

POO 

/ rft[7^u2 + c-%^(u ■ uf] < oo (7) 
Jo 

on the critical surface. This is possible only if the asymptotic condition 
lim u(t) = holds. Off the critical manifold u(t) diverges. Thus given 

r(0), r(0), the asymptotic condition singles out the unique r{0) on the critical 
surface. 

Inserting the asymptotic condition in (|1|), we see that — V0(r(t)) — » as 
t ^ oo, which implies in essence two distinct scenarios, (i) The particle is 
scattered into a region where F'^'^ = 0. Then u{t) has a limit as t — *• cxd and 
r{t) grows linearly, (ii) The motion is bounded. Then the particle comes to 
rest, lim u{t) = 0, at a point where the electrostatic force, — V0, vanishes. 

t— >oo 

Note that in general the condition — V0(roo) = does not determine the 
asymptotic position Too. E.g. a uniform magnetic field is confining even for 
0(r) = 0. 
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In analogy to (j^), our next task is to derive an effective second order 
equation for the motion on the critical surface. We follow the steps leading 
to and obtain 

+ (e/m)2F^"(x)F/(x)t;, + (e/mc)2F'^"(x)F/(x)t;<,t;,t;^} . (8) 

In principle one could compute also higher order terms. But they have the 
same magnitude as those contributions neglected already in the derivation 
of the Lorentz-Dirac equation. In addition (H) correctly describes the long 
time behavior as dominated by radiation reaction. Higher orders yield no 
qualitative change and, at best, make a minute correction of relative order 
10~^^ or even smaller in concrete examples. 

Eq. (H) appears in the second volume of the course in theoretical physics 



by Landau and Lifshitz |]TT|, who were guided by the insight that radiation 
reaction must have a small effect. One can only speculate why the Landau 
and Lifshitz equation (§) is apparently ignored in the literature. For sure, 
they do not discuss the structure of the flow with its critical manifold nor 
the relation to the asymptotic condition. 

There are several cases of interest where the solution to (|^) can still be 
handled analytically. The first one is a vanishing magnetic field and an 
electrostatic potential varying only along the 2-axis. Setting r = (0,?/,0), 
u = (0,?/,0), e0(r) = V{ii), the one-dimensional motion is governed by 

1 (^72/) = -V\y) - {e^/QT^^){l/m)V"{y)^y. (9) 

If V is convex, the energy is damped monotonically. The particular case of 
a quadratic potential is studied in the recent third edition of the textbook 
by Jackson |jl5| in the context of line breadth and level shift of a radiating 
oscillator. However, if V is periodic, say Viy) = Vocos{koy), then at the 
maxima the particle gains in energy from the near field, a process dominated 
by the energy loss at the minima. For long times the particle comes to 
rest. Also if V has a linear piece, then in this spatial interval the charge is 
accelerated without friction. 

An experimentally more accessible set-up is the motion in a uniform 
magnetic field (0,0,5). Then (||) simplifies to 

Y (m-fu) = eBu^ - {e'^ /Q7vc^){eB /m)^-f^u (10) 

(JjL 
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with u = {ui,U2,0) = u{cosip,smip,0), = {—U2,Ui,0). Setting a = 
e^/Gvrmc^ and the cyclotron frequency Uc = eB/m, one obtains du/dip = 
—aujcip, i.e. u{(p) = u{0) exp[—aucf]. For electrons = 8.8xlO~"^®i?[Gauss] 
which shows that radiative reaction is very small even at high fields. The 
charge spirals to its central rest point. Using ( |10|) the radius shrinks as 

r(t) = roe-"-^'V(l + (7 - 1)((1 - e-'^^'^') / 2)) . (11) 

In the ultrarelativistic regime, 7 3> 1, (|n|) simplifies to r{t) = ro(l + 
•yau'^t)'^, provided au^t <^ 1. To have some order of magnitude, in the 
case of an electron, au"^ = 1.6 x 10~^(-B [Gauss] )^/ sec and tq = 1.7 x 
10~^ (7/(5 [Gauss] )) m. Thus for B = 10'^ Gauss and an ultrarelativistic 
7 = 6 X 10^ the radius shrinks within 0.9 sec from its initial value tq = 10 cm 
to r = 1 fim, at which time the electron has made 2 x 10^^ revolutions. Only 
then the power law, (1 + crosses over to an exponential damping. 



Our third example is the motion of an electron in a Penning trap [12 
The electron is subject to a uniform magnetic field, as before, and in addition 
to the electrostatic quadrupol potential 

1 11 

ecpix, y,z) = - mul{--x^ - -y"^ + z^) . (12) 

A non-relativistic approximation suffices and in (H) we only keep terms to 
linear order. Then the in-plane and axial motion decouple and satisfy, with 
r = (a;,2/),u = (mi,U2), 

1 p2 1 1 

r = -a;,r + ^.u - {(^c - 3 ^> + 2^.^.r } , (13) 

e2 

z = -uj\z-- — —u:\z. (14) 

The in-plane motion can be solved easily. Without friction there are 
two modes with frequencies = \ [ooc ± — '^^zY^'^)- For the damping 
coefficients first order perturbation theory is sufficiently accurate with the 
result 

^+ = ' ^- = ' (1^) 

in agreement with a QED resonance computation [|T^ . Of course with some 
extra effort, one could handle also nonlinear and relativistic effects. We 
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emphasize that (JT^) is beyond the capabihty of Larmor's formula, which 
works only for a single mode. 

Note that uj- < 0;+ and therefore 7_ < which reflects that a loss of 
energy due to radiation lowers the potential energy and increases the orbit 
size. In practice lOcjl-n = 164 GHz and uJz/27i = 62 MHz. Inserting in (|15D 
leads to 7^^ = 8 x 10"^ sec and 7!^ = —3 x 10^^ sec. Experimentally, one 
observes that the magnetron motion (ti;_) is stable over weeks, whereas the 
cyclotron motion (a;+) decays to equilibrium within fractions of a second. 

To summarize, we have investigated the solution flow of the Lorentz-Dirac 
equation and discovered that in its structure it is identical to renormaliza- 
tion group flows in critical phenomena. The physical solutions are on the 
critical manifold and are governed there by an effective second order equa- 
tion. This equation is not plagued by the difficulties usually associated with 
the Lorentz-Dirac equation. In particular, the solutions to are stable 
and have the correct long-time behavior. Our examples show how radiation 
damping can be handled systematically and with ease. It would be of in- 
terest to have more stringent experimental tests. E.g. one could decrease 
the magnetic field in the Penning trap, so that the two modes mix better, 
and try to reach the resonance point u"^ = where the life-times should 
vanish according to (|l^). A further interesting possibility is to turn the mag- 
netic field out of the symmetry axis. Then all three modes mix resulting in 
damping coefficients which can be understood only on the basis of (^. 



*** 

I am most grateful to Fritz Rohrlich for instructive discussions and for the 
hint that the independently derived Eq. appeared in Landau and Lifshitz 
already a long time ago. I am indebted to Wolfgang Schleich for pointing 
out the Penning trap as an interesting application. I thank M. Rauscher for 
details on synchroton sources. 
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